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Abstract 



We introduce the notion of Multichannel Conformal Blocks relevant for the Operator Product 
Expansion for Null Polygon Wilson loops with more than six edges. As an apphcation of 

these, we decompose the one loop heptagon Wilson loop and predict the value of its two loop 
OPE discontinuities. At the functional level, the OPE discontinuities are roughly half of the 
full result. Using symbols they suffice to predict the full two loop result. We also present 
several new predictions for the heptagon result at any loop order. 



1 Introduction and Review 



Null Polygon Wilson loops in A/" = 4 SYM are extremely interesting objects. One of the 
reasons is because they yield all planar scattering amplitudes of the theory [1].^ In [2] an 
operator product expansion (OPE) for Null Polygon Wilson loops (NPWL) was proposed; it 
is the analogue of the usual operator product expansion for local operators in a conformal field 
theory. One of the main ingredients in the study of the latter are functions which package 
together the propagation of a conformal primary and all its descendants in a conformal field 
theory. These are known as conformal blocks [3, 4]. The main focus of this brief note is 
on the OPE for NPWL with more than six edges and in particular on their corresponding 
conformal blocks. They were first used in this setup in [5] in the re-derivation of the two 
loop hexagon Wilson loop [6, 7]. 



OPE review for the simplest kinematics 

Let us recall how to interpret simple known results for NPWL from the OPE picture. For 
more details see [2, 5, 8]. We start by discussing some kinematics before turning to the 
dynamics. The simplest NPWL are polygons which lie in a M^'^ plane [9]. The cusps of 
such polygons are given by light-cone coordinates (x^,a;,~). The conformal symmetry group 
acting on such polygons reduces to SL{2)^ x SL{2)^ generated by 



-1 ' 



dx] 



dxt' " ' dxt' ' ^ ' ' dxt 



with similar generator La acting on x^ . Conformal cross-ratios factorize into a product 

of left and right cross-ratios y^i,, = ^+!!^+ and v",,, = ^^J}^t ^^J. which are 

invariant under SL{2)-^- x SL{2)^ transformations acting on all cusps. Such transformations 
are generated by 

n n 

La = Y.Lf and I. = 5^I« (2) 

i=l i=l 

where 2n is the number of cusps. The first non-trivial polygon occurs when we have four 
x'l and four x~, that is for n = 4. This is the octagon represented in figure 1. It contains 
only two independent cross ratios which can be chosen as = 6"^"^ and x~ = 6^°". One 
configuration of xf leading to these cross ratios is depicted in figure la. The variables r 
and a have one further nice geometrical meaning as we now explain. Given a NPWL we 
can define a null reference square represented by the dashed lines in figure 1. A choice of 
square automatically splits the polygon cusps into two groups: the top and the bottom as 
represented in the figure. ^ There is a IR+ x ]R_ c SL{2)^ x SL{2)_ residual conformal 



^The usual bosonic loops, which are the ones considered in this note, lead to Maximally Helicity Violating 
amplitudes. A proper supersymmetrization of bosonic loop leads to all scattering amplitudes. 

^In the usual OPE for local operators we pick a conformal frame which consists of two points, conven- 
tionally chosen to be zero and infinity. The conformal generators used for the OPE (dimension and spin) are 
those that leave the conformal frame invariant. The reference square is the conformal frame in the Wilson 
loops picture. 
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Figure 1: The octagon NPWL in M}'^ kinematics, (a) The Euchdean picture where all points 
are spacelike or null separated. This is the kinematical regime considered in this paper, (b) 
A more physically intuitive Lorentzian picture. The results for the two pictures are related 
by analytic continuation. 

symmetry which leaves the reference square invariant. For our choice of points, the residual 
symmetry IR+ {^-) is generated by Lq (Lq). Consider the action of these symmetries on one 
of the groups, say the bottom one. That is, consider 

Co= Yl = -25- and Co = Yl 4*^ = 25- (3) 

iSbottom iSbottom 

We see that the action of the residual symmetry is trivial: it amounts to translations of r 
and a. An equivalent way of rephrasing this is the following. We start with any octagon, for 
example the one with t = a = 0. Then we generate a family of polygons by acting on the 
bottom cusps with the symmetries of the reference square. For the octagon in M^'^ we have 
two symmetries and two independent conformal-cross ratios. Hence we describe all possible 
octagons in this way. 

Having discussed the kinematics we move to the discussion of the Wilson loop expectation 
value. We will continue to base our discussion on the simplest example, the octagon in M^'^. 
Wilson loops with cusps have well understood UV divergences coming from virtual gluons 
exchange between neighboring edges. A regularization of these divergences breaks conformal 
invariance in a controlled way [10]. Since they are well understood it is simple to subtract 
them. What remains is a well defined conformal invariant function dubbed the remainder 
function. An alternative observable is constructed by taking a ratio of NPWLs such that the 
divergences cancel out and one remains with a conformal invariant finite function named r 
in [2]. It is related to to the remainder function in a trivial way [2, 8]. For the octagon, that 
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ratio of Wilson loops is 



/Tj/octagon\ /Ti/reference squareV 

^octagon ^ YofT ^ 1^- (4) 

^]/{/top hexagon^ ^yj/bottom hexagon^ ' ^ ^ 

see figure 2a. At one loop this quantity is given by the correlation function of the two 
disconnected Wilson loops as depicted in figure 2b. If two opposite cusps of the reference 
square are chosen to coincide with two cusps of the octagon (as in the figure), then all 
divergences exactly cancel.'^ We have [8] 

which is indeed a conformal invariant function of the two cross-ratios of the octagon. It is 
designed such that only excitations that propagate from the bottom to the top survive. 

Let us try to derive this expression from the OPE picture, that is without doing any 
Feynman diagram. We want to think of excitations being created at the bottom of the 
Wilson loop in figure lb, propagating in the middle region and being absorbed by the top 
part of the Wilson loop. In the propagating region the excitations interact with the left 
and right null edges of the reference square. Those two edges create a flux tube between 
them hence the excitations we have in mind are flux tube perturbations [2]. The flux tube 
excitations will have two important quantum numbers: energy (or twist) E conjugated to 
the r translation and momentum p conjugated to the a translation. To measure what is 
flowing we construct the Casimir of the S'L(2)+ conformal group, 

C = 2Cl- - (6) 

where 

i£bottom 

Note the following very important point. The generators Ca act on the bottom part of the 
polygon only. This is what we want to do since we want to measure what was created there 
and is propagating in the direction of the top region. In our parametrization we have 

The action of C on a primary excitation of twist E and its conformal descendants yields the 
casimir 2E{E — 2)} In the two dimensional kinematics, the excitations are just the primary 
F+_ and its descendents (D_)'^~^F_|__. These are the excitations which are generated by 
deformations of the loop in M^'^. The twist of the primary F+_ is = 2. Hence we should 



■^In other words, the total dual conformal anomaly [10] is zero for this quantity. 

^What we are doing mimics very closely the treatment of conformal blocks of Dolan and Osborn for 
correlation function of local operators [4] . 

^An important assumption is that the vacuum, represented by the flux, is SL{2) invariant; this is true 
at the first loop orders but eventually it breaks down [11, 8]. For the purpose of the current paper we do 
not need to worry about the higher loop breakdown. When studying higher loops it might be important to 
understand more precisely how this is broken. 
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Figure 2: Octagon ratio r. At one loop this ratio is computed by the disconnected correlator 
of two rectangle Wilson loops. In the OPE picture excitations are produced in the bottom 
rectangle and absorbed by the top rectangle; in between they propagate freely since at this 
loop order there are no interactions with the flux tube. That is, the gluon "doesn't know" 
about the reference square. 



find Cr'^^l^^°^{T,a) = which implies that r"^i^^p^{T, a) = ci((j)log(l + e"^"^) + 02(0"). The 
result should vanish in the OPE limit r — i- 00 and it should be cr -H- r symmetric (this is just 
parity, see figure 1). Hence 02(0") = and Ci{a) = clog(l + e~^'^). Up to the undetermined 
constant c related to the strength of the couphng constant, we just derived the octagon one 
loop result (5) from conformal symmetry. In this example, we see that 

T{r) ^ log(l + e--) = Y: ^^"'^ (9) 

k 

describes the propagation of the primary and all its descendents. In the usual OPE language 
we would rephrase this by saying that J^{t) is a SL{2)^ conformal block describing the 
propagation of a primary of twist E = 2 and all its descendents. It is fixed by the differential 
equation CJ^ij) = 2E{E — 2)J^(r) with the appropriate boundary conditions. 

What about higher loops? To think about what could change it is instructive to re-write 
(5) in the following inspiring form [2], 

rTLTi^, ^) = E / dpCi^^(P) e-^^--^^°'(^)^ (10) 

k 

where E^\p) = 2k is the free twist of F+„ and its descendents and the form factor C^\p) 
is interpreted as the probability amplitude for creating the excitation with energy at 
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the bottom and absorbing it at the top. It can be computed by Fourier transforming the 
one loop result. The momentum p is the quantum number conjugate to the translations 
in a which is a symmetry of the reference square. There are three types of higher loop 
corrections: (a) We can start having more than a single particle propagating,*^ (b) the form 
factors can get corrections and finally (c) the classical energy of the excitations will acquire 
an anomalous contribution from the interaction with the flux tube. The latter has quite a 
distinctive feature in perturbation theory: it appears in the exponent multiplying r. The 
expansion of the exponential in the 't Hooft coupling g"^, leads to a term linear in r, i.e. we 
have 

rTS:ir, cr) = Cf\^)e-'' +r ^ Df\a)e-'^ (11) 



from contributions (a) and (b) from contribution (c) 

The anomalous dimension of the excitations D^^F^ reads [8, 11] 

Ei\p) = 2/ [V^(l + ip/2) + ^{l - ip/2) - 2i,{l)] = g'^2{p) (12) 

In particular, it is independent of k. This is to be expected since these excitations are 
descendents, related to -F+_ by the SL{2) symmetry that is not broken at this loop order.' 
Hence we can easily compute the contribution linear in r of the two loop result,^ 



octagon 
' 2 loop 



-'-)|lincarm. = ^i^^P) 7^ (l^) (13) 



Similarly, at / loops the OPE expansion indicates that 



k 



+ ^^^E / dpCi'\p)^,{py~'e-^^-''- (14) 

This provides an infinite amount of predictions for the octagon Wilson loop at arbitrary loop 
order. At two loops this is enough to constrain the result completely [8, 5].^ At higher loops it 
provides important constraints but it does not seem to be enough, one needs to consider also 
other contributions coming from corrections to form factors, multi-particles, higher 

loop corrections to the anomalous dimensions etc. 



New ingredients for general kinematics 

Let us summarize the discussion so far and point out the new ingredients that come into 
play when considering more general polygons (they are the main focus of this short note) 

^We have gapped excitations on the infinite Une and therefore can talk about particles. 

^At two loops however, the SL(2) symmetry is broken and the two loop energy starts to depends on 

k. 

^This agrees neatly with the known result [12] as shown in [8]. 
^Which agrees with earlier predictions [12, 13]. 
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Figure 3: A five point function of local operators. It is the correlation function analog of 
the heptagon NPWL considered in this paper. For the five point function, the generalized 
conformal blocks describe the propagation of two primaries and their conformal descendants 
in two channels simultaneously. In the figure, the two primaries are parametrize by their 
dimension (A) and spin (/). 

• For the octagon in M}'^ we have 2 independent cross-ratios. For a polygon in M}'^ with 
n edges we have 3n — 15 independent cross-ratios. 

• For the octagon in M^'^ we could parametrize the 2 cross-ratios by the symmetries of 
the reference square. The main advantage of this parametrization is that it allows us 
to make use of the OPE physical picture in a straightforward way. For the hexagon 
in M^'^ we have three cross-ratios and this is still possible. For this we use the two 
symmetries that we used before plus an extra symmetry of the reference square which 
we did not make use so far: the S0{2) rotations in the two directions transverse to 
the plane where the square is. The square has no more symmetries. Hence we can not 
parametrize the family of all heptagons (or any other M^'^ polygons with more than 
six edges) in such simple geometrical way. In this paper we propose a parametrization 
of higher n polygons using more than one reference square. This will still allow us to 
make use of the OPE in a very efficient way. For example we will predict what the 
analogue of (13) and (14) is for the heptagon Wilson loop in M}'^. 

• For polygons in M^'^ the relevant primaries which we should consider at one and two 
loops are F+_ excitations. For polygons in M^'^ there are two infinite towers of primaries 
which we need to consider [2, 5] (they are labelled by the transverse S0{2) charge which 
played no role for the two dimensional kinematics). Hence, we will have to deal with 
an infinite tower of conformal blocks instead of the single block (9) considered above. 
For the hexagon polygon this was done in [5]. 

• Polygons in M}'^ with more than six edges are the analogue of correlation function for 
local operators with more than four points. The kind of conformal blocks that we need 
to introduce depends on more variables and capture the propagation in more than one 
channel simultaneously. For example, for the five point correlation function as well as 
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for the heptagon NPWL, these blocks describes the propagation of two primaries and 
their conformal descendants in two different channel at the same time (see figure 3). 
The computation of these blocks on the heptagon is our main result. It will allow us 
to decompose the one loop result in an OPE friendly way and therefore it will allow us 
to make an infinite number of predictions for the heptagon at any loop order. 

In this paper, we explain how to bootstrap polygons with more than six edges using the 
OPE program of [2]. The basic novelty of the present approach is the consideration of more 
than one OPE channel simultaneously. To illustrate the power of our approach, we apply 
it to the heptagon NPWL. We compute its multi-channel conformal blocks and use these 
to decompose the one loop result in an OPE friendly way. It will allow us to compute the 
heptagon OPE discontinuities and make an infinite number of predictions for the heptagon 
at any loop order. The paper is organized as follows. Section 2 is the bulk of this short note 
where we derive the heptagon multi-channel conformal blocks, decompose the one loop result 
and use it to compute the OPE discontinuity. In section 3 we discuss the generalization to 
NPWL with more edges. The appendix include some technical details for the heptagon. 

2 The heptagon Wilson loop 

In the two dimensional kinematics the first NPWL with non-trivial cross-ratios is the octagon. 
In four dimensional kinematics, the first non-trivial polygons is the hexagon. These are 
the analogues of the four point function in the context of correlation functions of local 
operators. This paper is about polygons with more edges which are the analogue of higher 
point correlation functions. 

In this section we focus on the heptagon in M}'^ which is the Wilson loop cousin of the 
five point correlation function of local operators. It contains all the new ingredients involved 
in the generalization from six edges to higher and hence serves well the purpose of explaining 
the general method. Our main result of this section are the decomposition of the heptagon 
one loop result in multi-channel conformal blocks presented in (21-23) and the prediction of 
its two loops OPE discontinuity presented in (26). The generalization to higher polygons is 
discussed in the next section. 

The first step in bootstrapping the heptagon is to choose an OPE channel and a corre- 
sponding conformal frame. That is, to choose a reference square. For the heptagon, all such 
choices are related by cyclic permutations. One then constructs the corresponding ratio of 
polygons (see figure 4. a) 



This is a finite conformal invariant observable. The function r is associated with a specific 
channel and therefore, contrary to the remainder function R, it is not a cyclic invariant 
function. The relation between the two is simply 



r = log 




;op hexagonV /T/t/bot pentagon 



;agon\ /Ti/red dotted square 



(15) 



'heptag( 



;on 



'top hexagon -'l-bot pentagc 



;on 



(16) 
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(a) (b) 

Figure 4: (a) Finite, conformal invariant ratio of Wilson loops using a single reference square, 
(b) Finite, conformal invariant ratio of Wilson loops using two reference squares. 



where r^^^^^ is the ratio (15) computed in a U{1) theory and dressed by the cusp anomalous 
dimension [2]. In other words, Tfj^^^ is the ratio (15) obtained by replacing the Wilson loop 
expectation values by the BDS result [15]. An advantage of r over R is that it is non zero 
already at one loop. At one loop it contains a single excitation propagating freely through 
the flux tube. 

The next step is to decompose r at one loop in terms of excitations propagating on the flux 
tube represented by the reference square. At one loop, the propagation of these excitations 
are organized in SL{2) conformal blocks. At this point we face a technical obstacle. The 
SL{2) conformal transformations form a three parameter group. The heptagon on the other 
hand has six independent conformal cross-ratios. In general all of them transform under 
the SL{2) symmetry which preserves the two null lines associated to the channel under 
consideration. That makes the decomposition very hard. 

Fortunately, there is a simple way around which - in addition - will also allow us to 
make many more predictions at higher loops. Instead of choosing a single channel for the 
OPE expansion, we choose two adjacent channels. We then study excitations that propagate 
in both channels. For the heptagon, any such choice is related by cyclic permutation. In 
figure 4b we plotted two reference squares, red (dotted) and blue (dashed), corresponding to 
the two adjacent channels. We will refer to these as the bottom (red dotted) and top (blue 
dashed) channels. Instead of the ratio r (15), we now consider the new finite conformal 



9 





-Li2 (1 - ui) - Li2 (1 - U2) - Li2 (1 - U4) 

-Li2 (1 - ^5) - Li2 (1 - U,) + Li2 (;^) 



1 



l0g(]ii^)l0g(l-U6) 



^ U2UQ — I 

- log (7/4) log (M5) + log (1 - U3) log (U1W5) 

-iogf --T:itr^ )i»g(i-"^"6) 



Figure 5: 
loops 



son 



At one loop r is given by the correlation function of two Wil- 
represented in the figure. The Ui are given by ui 7 = 



as 



^1,T^Zfi ■''l,3-''4,7 ■^l,5-''2,4 ■'•2,6-''3,5 •^3,7-^4,6 



2'"^ 2'^ , 2'*^ \\ only six of them are independent. 

"^1,5^4,7 ^1,5^2,6 J 



invariant ratio 

^yi^heptagon^ ^pj^ref pentagon^ 
-^Og ^yj^^Qp jjgxagon^ ^p^bot hexagon^ ' ('^''') 

That ratio is constructed such that it only contain excitations that propagate through both 
flux tubes represented by the red (dotted) and blue (dashed) squares. It is related to r and 
to the remainder function in a simple way that only involves lower point objects that we 
know already, 

^pp/bot hexagon^ ^pp/red dotted square^ 



r = r — log 



^pp/ref pentagon^ ^pp/top pentagon^ 



-^heptagon ~l~ ^[/(l) -^ref pentagon ~l~ -^top hexagon ~l~ -^bot hexagon • (-^^) 

Note that for r, there is no distinction between what we call top and bottom; it possesses a 
symmetry that exchanges the two parts of the polygon. At one loop, r is given by the cor- 
relation function of two Wilson loops as represented in figure 5. For the bootstrap program, 
we have to decompose the one loop result in excitations that propagate in both channels. 

Each of the channels is defined by a choice of two null lines and automatically defines a 
corresponding top and bottom group of cusps. Each channel comes with its three parameter 
group of SL{2) transformations which preserve the corresponding pair of null lines. These 
transformations are then taken to act on a group of cusps on one of its two sides. This is 
how we would measure what is flowing, see previous section and [4]. For example, we can 
choose to act with the S"/. (2) bottom transformations on the bottom cusps 2, 3, 4 (with respect 
to the red dotted square) and with the SL{2)top transformations on the top cusps 5, 6, 7 
(with respect to the blue dashed square), see figure 4. Since the two S'L(2)'s act on different 
cusps they clearly commute. In particular. 

The Casimirs of the two SL{2) 's commute. 

That is, we can decompose r into functions that are at the same time conformal blocks with 
respect to 5'L(2)bottom and 5'L(2)top. That fact is the main reason for studying r instead 
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of r. Without using these new blocks there would be no practical advantage; decomposing 
the new ratio r would have been as hard as decomposing the original ratio r. We will refer 
to these blocks as Multichannel blocks. The Multichannel blocks solve at the same time the 
differential equations arising from both Casimirs. With the appropriate boundary conditions, 
compatible with the OPE expansion, the Multichannel blocks are fixed uniquely in this way. 
Note that the Multichannel blocks, that are functions of the conformal cross-ratios, do not 
factorize into a conformal block of S'L(2)bottom times a conformal block of 5'L(2)top- That 
is because the conformal cross-ratios are constructed from ratios of distances of cusps on 
the bottom and top parts of the heptagon. Generically, all of them transform under both 
5L(2)'s.io 

To describe the Multichannel blocks we reconsider the two reference squares in greater 
detail. Each reference square gives us a three parameter family of polygons parametrized by 
the three symmetries of the square Rr^Ra^ S0{2)^. For the heptagon, we have two sets of 
these variables {r^, cTf,, 0b} and {rt, cxt, Together, they give us a basis of six independent 
cross ratios as needed for the heptagon. We see that a choice of successive OPE channels 
results in a natural basis for the conformal cross-ratios associated with these channels. A 
choice of these can be nicely represented using the Y functions of [16] as^^ 

_ (A5,Ai,A2,A3)(As,A7,Ai,A2) _ ^^—a^^—rf^ y- _ (A3,A4,A5,A2)(A4,A5,A6,Ai) _ ^^^^crt-n 

-^1.1 (A5,A6,Ai,A2)(Ar,Ai,A2,A3> "^2,1 (A3,A4,A5,A6>(A4,A5,Ai,A2> 



y- _ (A4,A5,Ai,A2)(A5,A6,A7,Ai) _ -2ti, y- _ (A4,A5,A2,A3)(A5,A6,Ai,A2) _ -2rt 

-^1.2 (A4,A5,A6,Ai)(A5,A7,Ai,A2) ^ ' "^2,2 (A4,A5,A6,A2>(A5,Ai,A2,A3> ^ 



Y _ (A4,A5,A6,A2)(A5,A6,A7,Ai) _ ^^^-cri,-Ti, Y = (A4,Al,A2,A3)(A5,A7,Ai,A2) _ (fn—at — Tt 

-^1.3 (A4,A5,A6,A7>(A5,A6,Ai,A2> ^ -^2,3 (A4,A5,Ai,A2>(A7,Ai,A2,A3> 

(19) 

The Multichannel conformal blocks are functions of the two SL(2) Casimirs Ci = — 1) 
as well as on the two momenta ki with respect to the a, directions. Here i = t,b stands for 
top or bottom. They are Appell hypergeometric functions of the second kind F2 [17] 

b t nj,=Orat=0 

r (- f + + /3b) r (- f + n, -f A) r (- f - f + + + /3, + A) 

r {n, + 213,) V {nt + 2A) 

where Tj = These Multichannel blocks are obtained by solving the corresponding two 
differential equation for the two reference squares as explained in [5]. They are normalized 
such that they are symmetric under the exchange of the top and bottom and are only 
functions of the Casimirs. That is, they are invariant under /3j — t- (1 — 

Note that the Multichannel blocks do not depend on 0^ and (pt in a direct way. Their 
only dependence on the corresponding charges is through [3b and /3t. To understand what 
values the /3's can take and what is their dependence on the two S0{2)^ charges mj, and rrit 
we have to discuss what type of primary excitations are propagating on the flux tube at one 
loop. These excitations were studied in detail in [11, 5]. For any 5*0(2) charge m, there are 



^"^In other words, the two differential equations, derived from the two Casimirs, do not admit solutions of 
factorized form. Instead, we have to solve both equations at the same time and the solution will entangle 
the cross-ratios in a non-trivial way. 

^"'^A choice of momentum twistors yielding these cross-ratios is given in the appendix 
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two types of excitations with /3 = m/2 and /3 = — m/2 corresponding to the two possible 
polarizations of the gluon. The type of excitation/gluon polarization does not depend on 
the reference square. Hence, an excitation with f3 = ±m/2 will have the same relation in 
any square. In our convention for the signs of (pb and (pt, it means that an excitation with 
Pb = ±^fe/2 will have Pt = l^^t/'^ correspondingly. Finally, note that the ratio r is parity 
invariant. Parity inverts the signs of (f)b and (f)f Therefore r depends on Q'nb<t>b+mt<i>t only 
through the combination cosh(mfc0fc + mt(j)t)- Given that content of excitations, parity and 
the symmetry of r under exchanging top and bottom, we are led to the general form of ru{i) 

mi,,mt=l 

where 

Am^^mtMM i4>b, 4>t) = Smt>l,mt>2 Cmt,mt-2,kt,kt COsh (m^^i + (m^ - 2)0fe) (22) 
+ Smt>2,mb>l Cmi,mt-2,kt,kt COsh {{mt - 2)0t + mfe^fe) 
+ Smt>2,mb>2 C2-mt,mt-2,kt,kt COsh {{mt - 2)(j)t + (2 - mb)0f,) 
+ Cmt~mb,kt,kb cosh (mt0t - mfc^fc) 

where 5x is equal to 1 if X is true and otherwise. In (21), we have used the symmetry 
of the SL{2) casimir C = (3{(3 — 1) under the exchange of /3 with 1 — /3 to represent all 
representation with /3 > 0. 

The one loop function for the heptagon is given in figure 5. It can be written indeed 
in the OPE form (21). Needless to say, this is a very non-trivial check of all the statements 
so far. To decompose ru(i) by brute force and hence read the form factors Cm^^m2MM is not 
an easy task simply because the function we want to compute is a complicated function of 
transcendentality degree two. However it turns out that we can construct some very useful 
box operators = (9|^ — d^^ and = — d^^. When acting on r, these operators simplify 
it dramatically.^^ For example the action of on ru(i) yields a rational function of the 
F-functions! This is highly non-trivial; naively we would expect lots of logs to survive and 
the expression to be a huge mess. Furthermore, further acting on this rational function with 
□t simplifies the result even more; once again, a priori, we would expect the exact opposite. 
We do not have a physical understanding of why these box operators are so remarkable, it 
would be very interesting to figure it out. At the end of the day all we need to do is to 
expand a (not so horrendous) rational function. To undo the action of the box operators we 
simply insert some propagator like terms. Then we check the decomposition numerically to 
be sure we did not loose any zero mode. The conclusion turns out to be remarkably simple, 

Cm,.^2MM = - _ ^ ^^^^ ^ rr^\ rj^ ^ • (23) 

^^This strategy was used in [5] as well to decompose the hexagon one loop result, see section 5.3. We see 
that the same kind of operators also simplify dramatically the heptagon result. This strongly indicates that 
there is something deep about them and that the Hexagon simplification was not a simple coincidence. 



12 



In (21) kh is integrated slightly below the real axis and kt slightly below the integration 
contour of ki,, that is kb EM. — iO and /c^ G M — 2i0. 

Given the one loop decomposition of r, it is straightforward to obtain a big part of the 
higher loop result. E.g., the two-loop remainder function can be split into two pieces as 

r2 loops = n [r2 loops] + [r2 loops] (24) 

The piece linear in r^, also called the OPE discontinuity, is obtained by dressing (21) by the 
one loop anomalous dimension 72si(fct) that is given by [8, 11] 



l2s{k) = -f2-2s{k) = 2g^' 



k k 



(25) 



Here, s is the conformal spin of the excitation. For /3 = \m\/2 in (22) we have s = 1 + /3 and 
for /3 = 1 + \m\/2 in (22), we have s = \m\/2. That is, [r2 loops]^^ is obtained from r[/(i) (21) 
by dressing A as 

Arni,'mt,ki,kMb,(t)t) 72+mt (/^i) '5mi>l,m6>2 Cmt,mb-2,fct,fc5 COsh (mt0t + (m^ - 2)0f,) (26) 

+ Irat ih) Smt>2,mt>l Cmt,mt-2M,kt COSh {{nit " 2)0i + mh4>b) 
Omt>2,mij>2 '~^2—mf,,mt — 2,kf,,kt 

cosh {{mt - 2)(j)t + (2 - mb)4>b) 

+ l2+mt {kt) Cmt~nib,kt,kb COsh (mj0t - mf,</)f,) 

More generally, at / loops there are / terms we can predict. ^"^ These are the terms that 
are proportional to Tj^^V^"^ with j = 1,2, ...,L They are obtained from r by dressing the 
excitations with 72st(fct)'~"' 7255(^6)"' as in (26). The corresponding terms in the remainder 
function are easy to read from (18). 

These are / new prediction at any loop order /. The next step involves a well posed 
mathematical problem. These predictions gives us part of a remainder function which we 
want to compute. We can repeat the computation using other OPE channels and hence 
obtain several constraints on the remainder function. The mathematical problem is to find a 
function with all the required symmetries and with the good OPE expansion in all possible 
channels. At two loops, if we assume that the answer has a symbol, one can show that the 
contributions computed above suffice to uniquely fix the function [2]. It would be interesting 
to see if these terms predicted here can fulfill the same role at higher loops. 

Finally note that our prediction expressed here as infinite sums should be equal to some 
function of transcendentality degree three. In [5] we succeed in resuming similar sums and 
got an explicit expression for the transcendental function. For the heptagon, we did not 
succeed in doing so. On the other hand, the sums yield a perfectly fine representation of the 
corresponding pieces of the answer; in particular these sums converge very fast (exponentially 
fast for large Xj). 



^■^It would be interesting to understand the Regge theory interpretation of these terms [19]. There seems 
to be a very interesting connection between the OPE expansion and the BFKL approach [18] as explained 
in [19]. 
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It is also possible to replace the decomposition mentioned above by the action of some 
convolution kernels and projection operators on the one loop result, see section 3.5 in [5]. 
The two approaches are equivalent, they are basically the Fourier transform of each other. 

Finally let us conclude with a discussion of the relation of our findings with those of 
[14]. In [14] the symbols of all two loop MHV amplitudes were computed. It would be 
very interesting to extract the seven point MHV amplitude from its symbol and compare its 
OPE discontinuity against our predictions. Alternatively, it would be great to compute the 
symbol of our sum, use it to construct the full symbol following [5], and compare it against 
the symbol derived in [14]. These two projects involve considerable work. 

We conclude with a third proposal which is considerably simpler to implement and which 
allows to "see" a great deal of the OPE expansion structure from the symbol of [14]. The 
proposal is to extract the symbol of the OPE discontinuity D2 from the symbol computed 
by Caron-Huot [14] and check that it satisfies the equation D^D_32 = 0. Here, D± = 
C — {d^ ± 2) projects out one of the two types of excitation, see [5] for more details. A 
more restricted check can be done using the multichannel structure. For the heptagon for 
example, consider the symbol of the combination r (18) and extract from it the symbol of 
the OPE discontinuity in either of the two channels. This is a simpler symbol D with three 
slots. Then, from the OPE decomposition and the relation (18) it follows that 

D^D^^'^B^ = DfD'^^^^ = Di'^Di^^Ds = D^^Z^f ©2 = 

where D^^ = Ci — d^. {d,j). ± 2) and the Casimir differential operators are written in the 
appendix. This third proposal has a clear advantage over the other two: it is straightforward 
to implement. Acting with differential operators on symbols is very simple, after all one of 
the definitions of symbols is through its differential, see e.g. [7]. The disadvantage of this 
check is that it does not probe the full OPE expansion. Namely, it probes the structure 
(21), (26) but it is insensitive to the form of the form factors. To probe the precise form of 
the form factors we have to work harder as described above. 

3 Generalization to more sides 

The techniques applied in the previous section to the heptagon generalize in a straightforward 
way to polygons with more edges. For an n sided polygon, one choses a sequence of n — 5 
successive channels as done in figure 4.b for the heptagon. That is, let the first channel 
on the "bottom" be associated with the two null edges {j,j + 3). Then the next channel 
be associated with the two null edges + 4) or (j — 1, j + 3) and so on. For a given 
bottom channel, there are 2"'~^ possible choices for the successive channels. For example, 
one choice is to have all the channels sharing the edge j. A natural symmetric choice 
is given by the Y function of [16]. That is the choice where the channels alternate as 
U,j + 3) (j, j + 4) -> (j + 1, j + 4) (j + 1, j + 5) . . . , see figure 6. a. Each of these 
channels comes with its SL{2) conformal group. The actions of all these S'L(2)'s commute. 

Next, for each channel one choose an associated reference square. A natural choice that 
generalize our choice for the heptagon is again in terms of the Y functions (see appendix D.2 
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Figure 6: (a) A particularly useful choice of successive channels convenient for describing 
NPWL using multichannel blocks. (b)Any null polygon Wilson loop can be decomposed into 
a sum of octagons, heptagons and hexagons, all sharing two null lines of the OPE expansion. 
The figure demonstrate such a decomposition of the top part of the loop (black) in these 
building blocks (blue thin lines). The two thick red lines represent the null lines of the OPE 
expansion. 

in [16]) 

Here, {tj, cTj, 0,} parametrize the three isometrics of the i'th square acting on the correspond- 
ing bottom group of cusps. The corresponding r for such sequence of squares is defined as 
in (17) using a reference n — 2 polygon. It is trivially related to the remainder function and 
lower n polygons as in (18). 

The Multichannel conformal blocks corresponding to the sequence of OPE channels are 
now functions of /cj, Tj} where i = 1, . . . , n— 5 parametrize the channel. They are obtained 
by solving n — 5 differential equations and are expected to be generalized hypergeometric 
functions of n — 5 variables. A decomposition of r in terms of the Multichannel block now 
allows for many predictions at higher loops. 

Given all the above, there is however a simpler way of decomposing the one loop result 
ru{i) for polygon with more then eight sides. It is based on the linearity of the one loop 
result and was first used in [8] to compute the two-loop remainder function for all polygons 
in M}'^ kinematics. The idea is to decompose any higher n polygon into a linear combination 
of octagons, heptagons and hexagons, all sharing two of their edges on the two null lines of 
the OPE expansion. That can be done as illustrated in figure 6.b. Due to the linearity on 
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the U{1) result, evaluating r[/(i) on the high n polygon is equal to the sum of the U{1) result 
on these smaller building blocks. In that way, a decomposition of the octagon, heptagon and 
hexagon automatically result in a decomposition of any higher sided polygon. 

Another way of getting the decomposition of ?"(7(i) for any polygon is by using the triangles 
of [5]. That is, r;7(i) is the one loop correlation function of two polygons as illustrated in 
figures 2.b and 5. As explained in [5], each of these two polygons can be decomposed as a 
sum of triangles. Therefore, the one loop correlator decomposed into a sum of correlators 
of pairs of triangles, one from the top and the other from the bottom. A decomposition 
of the contribution of such pair of triangle into excitations propagating in the flux tube 
automatically gives a decomposition of r^i) for any polygon. 

Both of these alternative approaches lead to a decomposition of the U{1) result in a single 
channel. They provide us with a restricted set of prediction for higher loops compared to 
the multichannel decomposition, see previous section. 

Finally let us end with a few words about the ground-breaking recent work done by S. 
Caron-Huot which appeared recently in [14]. In this paper the symbol of all MHV amplitudes 
at two loops was conjectured. We hope that the predictions presented in this paper for the 
remainder function will be useful to cross-check the results of [14] and in particular to help 
construct the Heptagon two-loop remainder function from its symbol. This was our main 
motivation for writing this short note. In the longer term we hope that the higher loop 
predictions presented in this paper will be useful in understanding the all loops result. 
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A Heptagon details 

In this appendix we present some details concerning the heptagon computation. It is of- 
ten useful to have a parametrization of the polygons in terms of momentum twistors. A 
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particularly simple such choice for the heptagon which leads to the cross ratios (19) is 
Ai = (1, 0, 0, 0) , A2 = (-1, 0, 0, 1) , A4 = (0, 1, -1, 1) , A5 = (0, 1, 0, 0) 

A3 = (-1, i^e-^t+at+^t ^ ^-2rt+24>tyl ^ _^rt-at-<Pt^ ^-n^at-^t ^ ^rt-at-4>t ^ ^ 

Ae = (0, 1, e^^^"^-'^" + 6-^*-'^''+'^''+"''-"*-"% 0) , A7 = (1, 0, e~'^''-"^+"^ e'^''"^-^^) 

Similarly, to compute the ratio (17) it is also useful to have a list of momentum twistors for 
the pentagon and for both hexagons, 

pentagon = {Ai, A2, A3+1/2, A5, A6+1/2} 

top hexagon = {Ai, A2, A3, A4, A5, A6+1/2} 

hot hexagon = {Ai, A2, A3+1/2, A5, Ae, A7} 

where 

A3+1/2 = (0, (1 + e<^*-^-'^*)-\ -1, 1) A6+1/2 = (0, 0, 1, 0) . 

There are two more polygons which appear in our discussion: the top and bottom reference 
squares. The mometum twistors of the bottom square can be parametrized as 

{Aleft, Atop, Aright, Abottom} bottom = {Al, (0, 0, 0, 1), A5, A6+I/2} (27) 

The Multichannel blocks are derived by acting with the two SL{2) Casimirs discussed in the 
text. Each SL{2) is the symmetry group which preserves two null lines. We have two SL{2) 
because we use two pairs of null lines, one pair for each of the reference squares in figure 4b. 
The action of each of the Casimirs is derived exactly as in [5], see section 5.2. Hence let us 
simply quote the final result. The differential equation for the bottom square reads 

(a - 4/3,(/3, - 1)) e-^^"'^"-*'^*'^* in, Tt) = (28) 

with 

C, = e-'^^ [drA, - drA, - d.A, + d^A, + (1 + e'^^X - 20^^, + + 2e'^A] 

The differential equation for the top channel is obtain from (28) with Cb — )■ Ct and Pb A- 
The two differential equations define the Multichannel block (20). In the two dimensional 
octagon case the Casimir operator was much simpler; it was simply given by (8). 

Concerning the decomposition of the 1 loop result, we find the remarkably simple result 

□ □ ^ 32^11^11^13121^22^23 p .29) 

* ' ""^'^ {YnYrsY22 + Yu (>2i>23 + (Via + >2i + 1) >22 {Yn + >23 + 1))) ^ ^ ^ 

where 

P = Fii (ri2 + l)((Fl3 + l) 122 + 1^21 (1^22 + 1))+ Vis (1^23 + 1^22 (>23 + l)) 

+ 1^12 (I2I (>22 + 1) (1^23 + 1) + (^-13 + 1) {Y23 + Y22 {Y23 + 1))) . 

It is probably possible to simply this expression even further using an overcomplete set of 
F-functions related to the basis used in this paper by the Y-system relations [16]. Expression 
(29) is much simpler to decompose than ^[/(i) given in figure 5. At the end, the action of the 
box operators is then simple to undo in Fourier by inserting some propagator like factors. In 
this way one derives (23). 
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